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We write a covariant transport equation for the phonon excitations of a relativistic superfluid 
valid at small temperatures. The hydrodynamical equations for this system are derived from the 
effective field theory associated to the superfluid phonons. We describe how to construct the kinetic 
theory for the phonon quasiparticles using a relativistic generalization of the analogue model of 
gravity developed by Unruh. This gravity analogy relies on the equivalence between the action 
of a phonon field moving in a superfiuid background with that of a boson propagating in a given 
curved space-time. Exploiting this analogy we obtain continuity equations for the phonon current, 
entropy and energy-momentum tensor in a covariant form, valid in any reference frame. Our aim 
is to shed light on some aspects of transport phenomena of relativistic superfiuidity. In particular, 
we are interested in the low temperature regime of the color fiavor locked phase, which is a color 
superconducting and superfiuid phase of high density QCD that may be realized in the core of 
neutron stars. 

PACS numbers: 47.75.-|-f, 97.60.Jd, 05.10.Gg,12.38.t,47.37.-|-q 

I. INTRODUCTION 

Superfluidity is a phenomenon that occurs in some systems at sufficiently low temperatures after the appearance 
of a quantum condensate [l|, . Landau gave an explanation of this effect and developed his famous two-fluid model 
considering that the system can be viewed as composed by a superfluid component, which describes the coherent 
motion of the condensate, and a normal fluid component, where dissipativc processes are allowed. The property of 
superfluidity follows from the existence of some elementary excitations with a linear dispersion law, the phonons. 
These phonons are the Goldstone modes which originate from the spontaneous breaking of a global symmetry caused 
by the appearance of the condensate. 

In the cold regime dissipative processes are mainly due to the collisions of phonons, which are part of the normal fluid 
component. The viscous hydrodynamics of a superfluid is involved because one can define more transport coefficients 
than in a normal fluid. However, for non-relativistic superfluids one can compute from the kinetic theory constructed 
by Khalatnikov [2] all the transport coefficients that enter into the two-fluid hydrodynamical equations. 

Relativistic generalizations of Landau's two-fluid model of superfluidity have been developed for applications to 
neutron stars physics. These generalizations are rather non-trivial, and alternative formulations of the hydrodynamics 
have been proposed in the literature [1, [3, d, 0, 111 ■ In the non-dissipative limit it is possible to show that all these 
approaches are equivalent. Unfortunately, very little is known about dissipative effects in relativistic superfluids [9l.[Tol. 

To the best of our knowledge, the whole set of dissipative terms which enter into the relativistic hydrodynamical 
equations has only been derived in one of the above mentioned approaches [§, [l^l ■ 

It is our aim to formulate the kinetic theory associated to the phonons of a relativistic superfluid, as these quasi- 
particles should dominate the transport phenomena at low temperatures. While the dynamics of the phonons in 
the frame where the superfluid component is at rest is easily formulated, this is in principle not so when done in an 
arbitrary frame. Relativistic boosts of the different variables complicate in an unnecessary way the problem. However, 
there is still the possibility to obtain transport equations in a fully covariant way . The approach exploits the gravity 
analogs which are naturally associated to relativistic hydrodynamical waves |12| , following ideas of Unruh |13| and 
others (see Ref.p^ for a review and references). 

Let us recall here that Unruh established long ago an analogy between supersonic flow in hydrodynamics and black 
hole physics [l^. Exploiting this hydrodynamics/gravity analogy [l^, Volovik pointed out that the kinetic equation 
of the phonons of a non-relativistic superfluid can be rewritten as that of a boson evolving in Unruh's acoustic metric 
[l6t . It is thus natural to expect that the phonons of a relativistic superfluid admit a similar treatment, as we show 
here. Let us mention that Popov also derived a similar transport equation for the superfluid phonons [Tj]- We will 
implement the same ideas in a somewhat different way. 

Relativistic superfluids may be realized in nature in the interior of neutron stars where the temperature is low 
and the energy scale of the particles is sufficiently high. In particular, in the inner crust of neutron stars the 
attractive interaction between neutrons can lead to the formation of a ECS condensate, the system then becoming 
superfluid. Moreover, if deconfincd quark matter is present in the core of neutron stars it will very likely be in a color 
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superconducting phase [31 . Quantum Chromodynamics (QCD) predicts that at asymptotically high densities quark 
matter is in the color flayor locked phase (CFL) [l9|. In this phase up, down and strange quarks pair forming a quark 
condensate that is antisymmetric in color and flayor indices. The order parameter breaks the baryonic number U{1)b 
symmetry spontaneously, and therefore CFL quark matter is a superfluid as well. 

The main part of the present paper is dedicated to the CFL superfluid. The starting point of our analysis is the 
low energy effective theory of the CFL phase which is obtained from the equation of state of quark matter [20| . 
Then we show how at asymptotic high density it is possible to derive from the microscopic physics, from QCD. the 
hydrodynamical behavior of the CFL superfluid. 

Our ultimate goal is having a well-defined formalism to derive all the transport coefficients associated to the 
CFL superfluid. Up to now only the shear viscosity and bulk viscosity associated to the phonon contribution of 
the CFL superfluid have been computed [U, as well as the contribution to the bulk viscosity due to kaons 
[2^, and an estimate of the thermal conductivity [23|. However there are other transport coefficients which are 
peculiar of superfluid systems that are still unknown. These are needed in order to derive the macroscopic behavior 
of a hypothetical compact star made of CFL quark matter, and single out possible signatures of quark matter in 
astrophysical scenarios. 

This paper is organized as follows. In Section |TT] we review the convective variational approach of rclativistic 
superfluid hydrodynamics. In Section IIIII we review the derivation of the hydrodynamical equations for the CFL 
superfluid phase of QCD and the low energy effective action for the supcrffiiid phonons. In Section HVl wc present the 
kinetic theory for the phonons of a rclativistic superfluid and derive the continuity equations for the phonon current, 
entropy and energy-momentum tensor. In Section |V] the transport equations are specialized to a local equilibrium 
state and the corresponding expressions of number density, entropy density and pressure are computed. We show that 
in the limit of small velocities these quantities agree with those of the phonons of a non-relativistic superfluid. We draw 
our conclusions in Section IVll In Appendix [XI we review an hydrodynamical approach to rclativistic superfluidity by 
Son. Wc will assume a physical Minkowski space-time metric (there is no real gravity field) and use the conventions 
^fj.i' _ j-^^ _i —1)^ with natural units c = ?i = fcs = 1. 



II. HYDRODYNAMICS OF RELATIVISTIC SUPERFLUIDS: CONVECTIVE VARIATIONAL 

APPROACH 

There are different formulations of the hydrodynamical equations governing a rclativistic superfluid [1, 0, H, [1, 0, H] ■ 
All of them were derived as rclativistic generalizations of the equations of Landau's two-fluid model of non-relativistic 
superfluid dynamics P, Q ■ Here we review the convective variational approach of Carter [J] . It has been shown in 
Rcf. that this approach is equivalent, at least in the non-dissipativc limit, to the potential variational approach 
developed by Lebedev and Khalatnikov Q. We also present in Appendix [X] the formulation derived by Son, which is 
equivalent to Carter's formulation as well. 

In the convective variational approach one defines a master thermodynamic function A, which is a scalar function 
of three independent scalar variables, namely n^np, s''sp and n^Sp, where is the entropy current vector and n'' is 
the (total) particle number vector. All the thermodynamic equations of the superfluid can be derived from A and by 
partial differentiation one finds 

dA = QpdsP + UpdnP , (1) 

that serves to define the so-called thermal momentum covector Qp and the vector Hp. 

By applying a Lengendre-type transformation, one can define the thermodynamics in terms of the pressure 

P = A - sPOp - nPfip , (2) 

which is a function of the conjugate variables = HpH'^, = ^pO'' and 8^ = 8''0p. Here /i is the chemical potential 
of the system, and O represents the temperature, which we will also denote by T . By partial differentiation one has 

dP ^ s^dBp + n^dup , (3) 

and then one can express the vectors and as a function of fi^ and as 

nP = FfiP + QOP , sP = QfiP + GO" , (4) 

where 

dP dP dP 

"-'W "^^a?- "^'^-^ 
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Finally the energy-momentum tensor of the system is given by 

TP" = nPn" + sPQ'" - PrfP" . (6) 

Although the notation used does not demonstrate the symmetry of the tensor explicitly, it really takes place 
The hydrodynamical equations of the superfluid are given by 

dpuP = , dpsP^O, (7) 

sPd^pQ^^ ^ , , (8) 

where the brackets denote index antisymmetrization. The irrotational condition equation on /ip tells that this vector 
can be expressed as the gradient of a scalar function, 

= ~dpV ■ (9) 

As it will be shown in the following Section, for the CFL superfluid it is possible to show that this scalar function is 
related to the phase of the diquark condensate [20l |. 

Using Eq. ([6]) and Eqs. ([7]) one obtains the energy-momentum conservation law 

dpTP" = . (10) 

From the expressions above it is clear that the hydrodynamical equations keep the form of conservation laws, exactly 
as in a normal fluid. 

In the zero temperature limit and for the ideal case of no interaction between the particles, the pressure is only a 
function of the chemical potential, and the entropy current vanishes. In this case, one defines the velocity vector 

aP 

vP = ^, (11) 
M 

such that it is properly normalized, vPvp = 1. In the T = limit the energy-momentum tensor takes the form of that 
of an ideal fluid and then Eq. ^ can be written as 

T"'3 = {nfi)v"v^^ - P?7"^ = iP + P) v^yf^ - Pt]"'^ , (12) 

where p is the energy density of the system, and we have used the zero temperature relation np, = p + P. 

Unfortunately in this formulation of the superfluid hydrodynamics it is not a priory obvious what is the physical 
meaning of all the hydrodynamical variables. The approach presented in this paper allows one to clarify their field- 
theoretical origin. 



III. THE CFL SUPERFLUID IN THE LIMIT OF VANISHING TEMPERATURE 



In this Section we study the relativistic superfluid system composed of color-flavor locked (CFL) quark matter. 
The CFL phase has long been known to be a superfluid: by picking a phase its order parameter breaks the baryonic 
number U{1)b symmetry spontaneously. At least in the limit of asymptotic high densities this is the energetically 
favored phase with respect to other less symmetric color superconducting phases and with respect to unpaired quark 
matter [l8|. 

In the regime of sufficiently high densities it is possible to derive the hydrodynamical equations from the microscopic 
physics. Indeed in Ref. [2^ Son showed how to obtain the effective action of the superfiuid phonons by integrating 
out from the QCD Lagrangian the heavy degrees of freedom corresponding to gluons, quarks and mesons. At very 
high densities, this is an allowed operation because the only low energy excitations are the phonons. The classical 
equations of motion derived from this effective action correspond to the zero temperature hydrodynamical equations 
that can be obtained within the convective approach, discussed in Section |TT1 This way of obtaining the Carter's 
equations allows one to identify the relation between the superfluid velocity and the gradient of the phase of the 
diquark condensate. Indeed, these two quantities turn out to be directly proportional. 

In order to see thermal effects in the low temperature limit in full generality, one has to consider the dynamics 
of the phonons moving in the superfluid background. This can also be derived from Son's approach (2^ . and we 
discuss in Section IIII Bl how to obtain the zero temperature effective action of the phonons in the presence of the 
superfluid background. As announced in the Introduction, the equation of motion of the phonon field in the superfluid 
background can be interpreted as that of a scalar fleld moving in an "acoustic" non-flat metric. This means that 
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we will write the equation of motion in a covariant form with respect to such a non-flat metric. This formulation is 
rather useful, as then one can infer the phonon effective action in any frame, and not necessarily in the superfluid rest 
frame. We will then use it in order to obtain the transport theory in covariant form. 

Let us stress here that we only consider the very low temperature regime, and that by increasing the temperature 
other quasiparticles can be excited in the C FL p hase. Although such excitations might be relevant for transport 
phenomena at moderately higher temperature [23|, for the sake of simplicity we will ignore them, as we assume a very 
low temperature. 



A. The effective Lagrangian of the superfluid phonons 

The effective field theory for the Goldstone boson originating from the spontaneous breaking of the U{1) b symm etry 
can be constructed from the equation of state (EoS) of normal quark matter, following the procedure of Ref. [2(| . 
Because the superfluid phonon is a Goldstone boson, and thus a light degree of freedom, its low energy effective 
field theory is obtained by integrating out the heavy modes from the QCD Lagrangian. This operation amounts to 
a minimization procedure with respect to the modulus of the order parameter of the symmetry. Exploiting other 
symmetries of the problem, one can entirely fix the form of the phonon effective action from the knowledge of the 
pressure of CFL quark matter. 

Calling if the phase of the condensate, and defining Dpip = dpip — fJ-qAp, with Ap = (1,0,0,0) and fiq the quark 
chemical potential, the low energy effective Lagrangian for tp (valid for energy scales much smaller than the the CFL 
gap A) is expressed as 

CMDpf] = P[{Dp^DP^y^^] . (13) 

In this expression P is a functional of the derivatives of the field tp that has the same form as the pressure of the 
system at zero temperature. Therefore in order to obtain the Lagrangian for tp is enough to specify P. At asymptotic 
large densities the EoS of CFL quark matter takes the form 

^(A^.) = i^/^^ (14) 

corresponding to a system of three flavors of massless quarks. The validity of such an expression relies on the fact 
that at very high /ig the coupling constant is small g{fJ-q) <C 1 and the effects of interactions and the effects of Cooper 
pairing are subleading. Therefore they are neglected in Eq. (|14p . One assumes that the quark masses give a subleading 
effect as well, because niq <C i^q. Effects due to the strange quark mass are considered in Ref. [22|. In the end one 
obtains from Eq. that the low energy effective Lagrangian for the CFL phase in the limit of asymptotic densities 
is given by 

JOes^^[{^o^-^^,r~i^^^?T■ (15) 

There is an interesting interpretation of the equations of motion associated to ip. Since the Lagrangian in Eq. (I15p 
does not explicitly depend on the field (p, but only on its derivatives, the corresponding classical equation of motion 
takes the form of a conservation law. Then one can view it as the hydrodynamical conservation law of a current 
representing baryon number flow, 

d^inov") = , (16) 

where 



dP 

no = -r- 



\fl' (17) 



is interpreted as the baryon density [S^, where fi = {DpipDPtpY^^ and where 

Dpip 
Vp = — , 



(18) 



is the superfluid velocity with (p the solution of the classical equation of motion. Notice that the superfiuid velocity 
is properly normalized, that is, VpV^ = 1, and that in this definition it is taken into account that Lorentz symmetry 
is explicitly broken by /iq. 
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The energy-momentum tensor can be written in terms of the velocity defined in Eq. (|18p and Noether's energy- 
density po, 

= MvPy'^ - gP'^Po = (po + Po)vfv^ - v'"'Po , (19) 

where we have written po + Pq = nojl, with Pq the quark pressure evaluated at /2. The energy-momentum tensor is 
conserved 

dX" = , (20) 

and traceless Tg^ = 0. 

One can immediately check that Eqs. (fT6|) and ([20]) are equivalent to the hydrodynamical equations of a relativistic 
superfluid at zero temperature obtained in the convective variational approach. The definition of the superfluid 
velocity differs in the zero component, but this can be matched by the change ip ^ ip + p,xo- 

B. Phonons moving in the superfluid background 



From the expression of the Lagrangian in Eq. pSj) it is possible to derive the effective field theory of the phonons 
moving in the background of the superfluid [22]. In the CFL phase the superfluid phonon is the Goldstone boson 
associated to the breaking of the U{1)b symmetry and it can be introduced as the phase of the diquark condensate. 
However, according to Eq. fTSl) the gradient of the phase of the condensate deflnes the superfluid velocity as well. Then 
it should be possible to decompose the phase of the condensate in two fields, the first describing the hydrodynamical 
variable, the second describing the quantum fluctuations associated to the phonons. Therefore in order to find the 
phonon dispersion relation in a moving superfluid we will consider the quantum fluctuations around the classical 
solution ip of the equations of motion associated to the Lagrangian (|15p . thus we write 



^(a-) - (^(.t) + 0(x) . (21) 

This splitting implies a separation of scales - the background field (^(x) is associated to the long-distance and long-time 
scales, while the fluctuation (f>{x) is associated to rapid and small scale variations. The gradient of (fix) is proportional 
to the hydrodynamical velocity, according to Eq. (jl8[) . while the fluctuation cj) is identified with the superfluid phonon. 
From the low energy effective action of the system 

S[^] = J d^xC,s[d^], (22) 

we deduce the effective action for the phonon field expanding around the stationary point corresponding to the classical 
solution 



4 d Cc« 

X ■ 



d{d^ip)d{d^ip) 

and considering the expression of the Lagrangian in Eq. (fT5|) one has that 



d^(l)d,<P+--- , (23) 



d{d^Lp)d{d^tf) 



"° + f ^ - 1 ) ^'^""^ ■ (24) 



Here Cs = l/'s/S is the speed of sound in CFL quark matter. 

The second term on the right hand side of Eg. ((23)) represents the action of the linearized fluctuation - here the 
superfluid phonon - and can be written as the action of a boson moving in a non-trivial gravity background [T^ 

S[4>] = (25) 

where we have defined the metric tensor 

= f]""" + - l) (26) 



and the determinant Q = l/dct|5''^|. 

In writing the effective Lagrangian in Eq. (|25p we have rescaled the phonon field according to 



(27) 



to get a dimensional scalar field, and thus a dimensionless metric in natural units. In doing the rescaling, we have 
assumed that on the scale of variations of 0, both /2 and no can be assumed to be constant. This is justified by the 
separation of scales which is implicit in the splitting of the bosonic field ip done in Eq. (I2ip . 
The normalized field obeys the classical equation of motion 

dJ^G'^'d.A =0. (28) 



We then see that from the effective Lagrangian in Eq. (jTS]) one can derive in Eq. (pS)) the so-called sonic or acoustic 
metric tensor f/'^^ [3] that depends on the speed of sound Cs and on the superfluid velocity v^^ . Note that since 
the superfluid velocity v'^ is normalized to 1, in general the metric tensor Q^'^ corresponds to a strong gravitational 
field. The limit of weak gravitational field is achieved for Cs — + 1 , and only in this (unphysical) limit one can expand 
the gravitational metric tensor around the fiat metric 77^^ . 

The dispersion law of phonons in the superfiuid medium are obtained by solving the equation of motion for a 
massless particle moving in curved space time, that is 

= , (29) 

where = {E, — p). In the superfluid rest frame, i.e. where v'^ — (1, 0), the dispersion law simplifies to the form 

E^csp, (30) 

which expresses the fact that phonons propagate at the speed of sound. In a different frame the dispersion law takes a 
rather complicated form. In order to express such a relation in a compact way it is convenient to define the following 
Lorentz invariants 

e = p^v^ , = - {rr - W)p^p, , (31) 

which arc the energy and the square of the 3-momcnta in the local rest frame of the superfluid. It is easy to check 
from Ea. ip^l that in any frame the relation e = CsTt holds. 



IV. KINETIC THEORY FOR THE SUPERFLUID PHONONS 



In order to study the contribution of the superfluid phonons to the hydrodynamics one has to consider a thermal 
bath of these quasiparticles, and compute their contribution to the various thermodynamical quantities. This analysis 
could be done in thermal field theory, using the action given in Eq. (j25|) . Because we will be mainly interested in 
studying kinetic phenomena, and ultimately in computing transport coefficients, we find more convenient to develop a 
transport theory for the cold CFL superfiuid. While it should be fully equivalent to the thermal field theory approach, 
it turns out that the computation of transport coefficients is much simpler when done with the use of kinetic theory, 
as it is the case for other field theories (see, e.g. psl]). 

Transport equations can only be formulated when there are well-defined quasiparticles in the system, that is, 
thermal excitations that are sufficiently long lived. In the CFL superfiuid a direct computation of the damping rate 
for the superfiuid phonons shows that it is suppressed by {T/^q)^ with respect to their typical energy [2l[. One can 
thus guarantee that it is possible to write a transport equation for these quasiparticles, as quark matter is in the CFL 
phase only when T <C . 



A. Quasiparticle picture 

As described in the previous Section, the dispersion law of the phonon moving in the superfiuid background can 
be seen as that of a boson propagating on a gravity background characterized by the acoustic metric. In a classical 
approximation one can view the phonons as massless quasiparticles that follow the paths determined by the geodesies 
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of the acoustic metric, defined by = Q^^^dz'^dx'^ . Moreover, it is possible to define a quasiparticle Hamiltonian, H, 
in sucli a way that solving the Hamiltonian equations one obtains the geodesies of the problem. Thus, we introduce 

n = \q^''p^p. = Iq^up^p" = ip^P^ , (32) 

where p^ — Gfj.iyP^ and where G^iiy is the inverse of the metric G'^'^ , obtained by solving the equation G^^Gvp = 5'^- For 
the metric tensor defined in Eq. (|26p the inverse is given by 

Gij.u = tltiu + {cl - l) WpWi. , (33) 

where ~ "n^j^^jV^ . 

The Hamiltonian equations of motion are given by 

however, since we will formulate the transport approach using the variables {x^ ^p^), rather than {x^,pp), it is conve- 
nient to find the equation of motion for pP. After simple algebraic manipulations one finds 

^ = -G'^iG.fj.. - ^tJ./3,p)pV = -r>V , (35) 

where F^^ denote the Christoffel symbols associated to the metric Gfj.i^, 

Ta/J ^ {Qt^P.c + Gf,a,P - • (36) 

With these Christoffel symbols one can define the covariant derivatives for a generic controvariant vector A'' as 

Af, = Af^, + r(;„A" , (37) 

where A^^, = d^A^^ . While for a covariant vector one has 

A^-,y = A^^i, — V^j^Aa ■ (38) 



B. Covariant transport equation for the superfluid phonons 

We will now study the Liouville equation governing the evolution of the phonon distribution function f{x,p). Here 
it is understood that f{x, p)p^ny,dTi dV is the number of particles whose world lines intersect the hypersurface element 
UfidTi around x, having four- momenta in the range {p,p + dp), with a four light-like vector that we will take as w^, 
and the momentum measure dV is chosen in such a way that it is coordinate invariant (26j 



dV = V^2H{p)SiG^.p'Y)j0^ , (39) 



where H{p) = 1 if p is future oriented for an observer moving with velocity w^, and otherwise. 

The presence of the superfluid background will be manifest in the fact that the Liouville equation has to be written 
in a covariant form with respect to the metric G'^'^ ■ Indeed, the evolution of the distribution function is governed by 
the equation 

where C[f] is the collision term and upon using Eqs. (|34I35|) the Liouville equation can be rewritten as 
that is the general relativistic version of the Boltzmann equation (26j . 
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Knowing the distribution function one can construct the current, energy-momentum tensor and entropy associated 
to the phonons, which are given by the expressions 



T' 



a0 
ph 



'ph 



[flnf-il + f)lnil + f)]dr 



(42) 
(43) 
(44) 



Notice that in Section UlIB I we expUcitly made a spUtting of fields, separating the low energy modes from the high 
energy modes, the last corresponding to the phonons. However, when writing the transport equation for the phonons 
one should keep in mind that the collective and thermal motion of the phonon fluid can occur at long scales. For this 
reason one has to keep the Christoffel symbols in the transport equation, which contain derivatives of the superfluid 
velocity. As it will be shown in Section |Vl the equilibrium distribution function of the phonons depends on Q^^ , 
meaning that all the above quantities evaluated at equilibrium vary on scales of the order of the scale of variation of 
the acoustic metric. 

From the expressions of the current, energy-momentum tensor and entropy one can derive the continuity equations 
obeyed by these quantities. The analysis is quite simplified if one takes into account the following property [26| 



pMi ...pl^ip^f(x,p)dV 



p"' ■■■p^"L[f]dV. 



(45) 



In the case of vanishing collision term, considering different values of I in this equation enables to obtain the various 
covariant conservation laws. 

Taking I = one has that the covariant continuity equation for the phonon number is given by 



(<h); 



p''f{x,p)dV 



L[f]dV = / C[f]dV 



which one can write explicitly as 



9.<h + r>;;i,= / c[f]dv, 



(46) 



(47) 



and expresses the fact that collisional processes can violate phonon number. Moreover one has to consider interactions 
of the phonons with the superfluid background. This is reflected in the second term in the l.h.s of the above equation 
where a Christoffel symbol appears signaling that propagation of phonons is taking place in "curved" space-time. 
Since we can write 



1^ 



(48) 



it follows that in the coUisionless case, the variation of the number of phonons stems only from a possible space-time 
dependence of the speed of sound. Therefore in the coUisionless case the covariant conservation law of the phonon 
current can be written in a compact form as 



("ph) 



1 



1 



(49) 



Notice that in the absence of collisions the quantity that is conserved with respect to the flat metric, i.e. with respect 
to the derivative da, is the current 

n^h = V^«ph- (50) 
This is the quantity that should be identified with the conserved phonon current, and obeys the continuity equation 

Sa^ph = , (51) 



in the absence of collisions. 
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From Eq. (j45p with 1 = 1, one obtains the covariant continuity equation for the energy-momentum tensor 

iT^^);.- j P^C[f]dV = Q, (52) 

where we have assumed that colhsions conserve both energy and momentum. The equation above can be written 
more exphcitly as 

(7^p7);. = iT^i:),^ + TL^pT + K.T;: = , (53) 

that can be rewritten as follows, 



[v^T^:) + v^T^i^n, = . (54) 

As in the case of the phonon current, it is useful to redefine the energy- momentum tensor as 

f^"' = ^/^T^"' , (55) 
and in this way the four-momentum associated with it is given by 

= f d^xV^T'^" = f (i3xf°^ (56) 



Notice that defining the energy- momentum tensor of the phonons as T^", we have that the total conserved tensor of 
the system is given by 

= T^'' + f^^ , (57) 

where T^^ is the component of the energy-momentum associated to the superfluid background. 

If we assume that the collisional integral is zero, one can expect to have also the covariant conservation of entropy, 
i.e. 

(s^h);a = -^dc (V^sl,) = -^d^ (S^h) = , (58) 



where Spj^ = \/ —Qs'^yi- I^- the presence of collisions this relation is not satisfied, and then there is entropy generation. 

Let us remark on the meaning of the continuity equations of phonons that we have derived Since the phonon 
fluid moves on the top of the superfluid, even assuming that the collision term vanishes we still have the possibility 
that processes of interaction between the phonon fluid and the superfluid will lead to the exchange of momentum 
and energy among them. Indeed, even if energy and momentum of the whole system composed by the superfluid 
condensate and by the phonons is conserved, there can be energy and momentum exchanges between these two 
subsystems. This is also what happens when considering matter propagating in a gravity field, as also in that case 
the energy-momentum tensor associated to the "matter" field is not strictly conserved. 

The possibility of exchanges of energy and momentum mentioned above is reflected in the continuity Eq. (|54p that 
we rewrite in the compact form 



d.T^x.. = ^T^d^g,,, (59) 



using the definition in Eq. ((55|) . On the right hand side of this conservation law there is a term that couples the energy- 
momentum tensor of the phonons with the derivative of the gravitational field, expressing in a covariant form the 
interaction between the phonons and the underlying superfluid. As a result, energy and momentum can be exchanged 
between the phonons and the background if the speed of sound is space-time dependent and/or if the superfluid 
velocity is space-time dependent. This fact can be understood also within the gravitational analogy. Indeed, the 
energy-momentum tensor of any matter field, here the phonon field, couples with the gravitational field G^^^, inducing 
exchange of energy and momentum between the gravitational and matter fields. 

Let us finish this Section by stressing that the interesting point of the transport theory presented here is that the use 
of the gravitational analogy automatically allows to write the equations governing the variations of phonon currents 
and energy-momentum tensor in terms of covariant equations with respect to the acoustic metric. 



10 



V. THERMAL EQUILIBRIUM OF THE SUPERFLUID PHONONS 



The covariant continuity equations for the phonons derived in the previous Section take a very simple form when 
the macroscopic variables associated to the superfluid are homogeneous or constant. Indeed in this case the ChristofFel 
symbols vanish, i.e. F^^^, = 0, and covariant derivatives are equal to ordinary derivatives. Thus, the kinetic equation 
for the phonons are the same that one obtains for ordinary particles evolving in a flat metric. In this case one can 
find the equilibrium distribution function for the phonons in the exact same way as for ordinary particles. 

One can as well find global equilibrium solutions for stationary situations (no time dependence). Global, or coUi- 
sional equilibrium, distribution functions are given by 



1 

exp {pi^f3i_, 



fcqix,p) = ,^^„^ , 7 , (60) 



where (3^ is a vector that can be derived as follows. Since i[/cq] = one has that obeys the equation 

(/3A,p-/3anp)pV = 0, (61) 

that can be rewritten as 

Px;p + Pp.X = . (62) 

Assuming that (3^ is time independent, it is possible to check that one possible solution to the above equations is 
given by 

(l/r,0) , (63) 

where T is a constant that we identify with the temperature. 

There are other local equilibrium configurations that are found by demanding no local entropy generation, that 
is sf^ = and with solutions given by the coUisional invariants of the theory. In this case we can expect that the 
equilibrium distribution function takes the form 

/eq(x,p) = } , (64) 

exp (p^/3^) - 1 

where 

= lii^xyix). (65) 
In the local equilibrium configuration the phonon energy-momentum tensor has to take the form 



= Hu^'u'' - KG"'' , (66) 



where H and K are two scalars. The reason for this form is that T^j^ is a tensor and because of the definition in 



Eq. (|43p it can be built by the tensors u'^u'^ and C/^" only. We will fix the scalars H and K as follows. When both 
the superfluid and phonon fluid are at rest, i.e. = = (1, 0), the pure temporal component of this tensor should 
give the phonon energy density at rest, while the spatial components are diagonal and proportional to the pressure 

rp"h" = Pph, T%^P^S^= . (67) 

Thus one can identify 

H = pph + Pph/c^ , K = Ppi,, (68) 
and the phonon energy-momentum tensor is then expressed as 

TpT = (Pph + Pph/cDu'^W' - Pphr;'''' - Pph (^^ - l) v^v'' . (69) 

From the definition in Eq. (|43p one sees that the energy-momentum tensor is tracelcss, i.e. Gp,i,T^^ — 0. Upon 
substituting the expression of T^^ that we have obtained above, one has that the traceless condition leads to the 
equation 

^ Pph f 4ci-g,.u^u'^ \ ,„^. 
'''' ^ ^ I g,.u.u^ ) • ^''^ 
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Moreover, using the definition of the energy-momentum tensor in Eq. ((55)) and the condition Q^uT!^^ = 0, one finds 

^ph = 4PphU^^" - PphG^" , (71) 

where we have defined 

Pph = V^Pph and (72) 



We also define — Q^i.vu'^ and therefore u^w^ ~ 1, meaning that is a four- vector normafized to 1 with respect to 
the metric 

It is important to note that the phonon contribution to the energy-momentum tensor keeps the form that is expected 
for the normal fluid component in a superfluid. Indeed, Eq. (|69p has the tensorial structure which is expected in these 
systems according with the Poisson bracket approach to the hydrodynamics reviewed in the AppendixIXl (see Eq. (|A2[) ). 

In a similar way one finds that the expressions of the current and of the entropy current are of the form 

"ph "phw^ , s'^i, = SphM^ , (73) 

where riph and Sph are scalars. 

Upon employing the definitions above in the Eqs. (|42)) . (|43)) and (|44p . with the local equilibrium distribution function 
()64p one can compute the phonon current, entropy current and pressure. Then the energy density can be deduced 
from the pressure, according to Eq. ([70]) . The best strategy for evaluating the integrals in Eqs. (j42|) . ([43]) and ([44]) 
is to choose the frame where the superfluid is at rest v^^ = (1,0). Then, after deflning u^^ = 7(1, u), where 7 is the 
Lorentz factor, we obtain for |u| < Cg 

^ _ T^Ci^) 1 ^ _ TV 1 p _ T^TT^ 1 

" ^2 c3^4(i_ ' ~ 45 ^374(1 - ^)2 ' - 90 c37^(l - ' 

After realizing that in the superfluid rest frame 

,2 



u 



g^.u^^u^ ^ciY{l~—) , (75) 
one can express the same quantities in any arbitrary frame writing 

"p'^""^ ^ {g^^u^^u'^y 'p'^"^ ^ 45 {g^^^u^^u-Y' ^^'^""^ ^ 90 {g^.w^u-y ^"'> 

From these expressions one can recover the non-relativistic results of Ref. 0. Indeed, taking the limits — > (1, ujvfl) 
and (IjVTVij), and using the same symbol for the non-relativistic speed of sound, one obtains that 

g^.u^'u" ^cl- {UNR - vnrY , (77) 

and the expressions above turn into the relations 

. _ T^CjS) 1 . _ TV 1 -3 _ T^TT^ 1 

"Ph - ^2 c3(l-^)2 ' ""P^ " 45 c3(l-^)2 ' ' 90 cf(l-^)2 ' ^ ^ 

where w = UAr_R — vnr is the so-called counterflow velocity 0- 

From the Eas. (|75|) one can also verify that the thermodynamical relation 

- /3^Pph + /3.Tp7 , (79) 



s 



ph ~ H ^ pii ^ Mi'-^ph 

is satisfled. In this relation indices are lowered/raised with the acoustic metric. 

A similar expression for the phonon contribution to the pressure in a cold relativistic superfluid was found by Carter 
and Langlois The main difference of our treatment with that of Ref. is that we assume that the velocity of 
the phonon fluid is relativistic. Carter and Langlois assumed a non-relativistic velocity for the phonon fluid, and 
accordingly, they obtained a non-relativistic relation among the thermodynamical variables instead of the relativistic 
one reported in Eq. (|75|) . 
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For completeness we write the conservation law of the energy- momentum tensor in Eq. (|59p in a more explicit form. 
Using the definition in Eq. (fTTj) we can rewrite the left hand side of Eq. ([55)1 as 

^M^'p';. = 4[(a^Pph)u'^?i'^ + Pph(a^^^)w'' + Pph(a^^i")w1 - a.Pph , (80) 

while the right hand side of Eq. ((59|) turns out to be 

^TpX d^Gp^ - iPph(4K^K'' " Q^d.Qp^. ■ (81) 
Projecting these equations along we obtain that 

4(a^Pii^) - ti^9^Pph = ^{^u^uP - gf^nu^d^Gp^ . (82) 
Projecting in the directions orthogonal to u'^ with {u°'u^ — Q"") one has 

- 4Ppi,g^''u''d^K - (*i"ii" - G'^nd.Pph - ^i^u^u" - ^^'')(u"u'^ - e"'')a,ep^ . (83) 

In the limit where the metric C/^j/ can be approximated with the flat metric 77^1/ , one has that — > and the 
equations above reproduce the hydrodynamical equations of a ultrarclativistic system. 

VI. DISCUSSION 

The properties of superfluids at low but non-vanishing temperature are influenced by the presence of phonons, the 
massless Goldstone bosons originating from the spontaneous breaking of a global symmetry. In case these quasipar- 
ticles are the only low energy degrees of freedom they will give the dominant contribution to the transport properties 
of the system in this cold regime. 

The transport theory of a non-relativistic superfluid was developed long ago, see e.g. Ref. 0. In the non-relativistic 
case the transport equations for the phonons can be expressed as those of a boson evolving in Unruh acoustic 
metric [l6j . In the present paper we have extended this approach to relativistic superfluids. We have shown that 
in a relativistic superfluid the same gravity analogy can be employed for deriving the transport equations replacing 
Unruh's metric with a generalized relativistic acoustic metric. Similar results have been reported in Ref. [l7j , although 
implemented in a different manner. The advantage of this formulation is twofold. It allows us to obtain the continuity 
equations and to express the phonon hydrodynamical variables in a covariant form, valid in an arbitrary reference 
frame. Moreover, it clarifies the physical meaning of the variables used in the convective variational approach of the 
hydrodynamics. 

The computation of transport coefficients can be implemented employing the framework presented in the present 
paper considering the Boltzmann equation (j4ip with a non-vanishing collision term. The relevant collision terms can 
be derived from the microscopic physics and depend on the particular system considered. For the CFL superfluid in 
the high density limit the collision term can be evaluated considering various scattering processes among superfluid 
phonons whose vertices can be read from the effective Lagrangian in Eq. (jl5|) . The leading processes are binary 



collisions, collinear splitting and joining processes [21[. The phonon contribution to the shear viscosity, and the bulk 
viscosity associated to the normal fluid component have been computed following this procedure in Refs. [U, [2^ . 
These computations were done using kinetic theory in the superfluid rest frame, where one takes w'' = (1,0), and 
assuming homogeneity in the superfluid flow. This corresponds to consider in our equations vanishing Christoffel 
symbols, P^j^ = 0. The transport equation that we have presented will allow us to compute the viscosity coefficients 
that involve also gradients in the superfluid flow. The results of these computations will be reported soon. 

Although we have focused our attention on the CFL superfluid the techniques developed here could be employed 
for different relativistic superfluids. In particular, it would be interesting to compute the phonon contribution to the 
different viscosity coefficients in the neutron superfluid that is supposed to be realized in the interior of neutron stars. 
The evaluation of all the bulk viscosity coefficients arising from non-equilibrium beta processes has been recently 
reported in Ref. [Toj . It would be interesting to evaluate the phonon contribution to these viscosity coefficients. One 
should also consider that at densities achievable in the core of neutron stars the energetically favored phase might 
not be the CFL phase [l^- Among the various possible phases that can be realized one interesting possibility is 
that quark matter is in the Crystalline Color Suprconducting phase [27| . where the gap parameter is periodically 
modulated in space and therefore spontaneously breaks translational invariance. In this case the low energy degrees 
of freedom do not consist only of the phonon excitation associated to the spontaneous breaking of U{1)b- Different 
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phonon excitations arising from the spontaneous brea king of space symmetries and describing the oscillation of the 
crystalline structure are part of the spectrum as well |28| . Moreover, gapless fermionic excitations are present and 
this makes the analysis of the low energy properties of this phase challenging. In such a situation the transport theory 
derived in the present paper should be extended to include all the remaining low energy modes. 
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APPENDIX A: POISSON BRACKET APPROACH TO THE HYDRODYNAMICS OF RELATIVISTIC 

SUPERFLUIDS 

The hydrodynamical equations of a relativistic superfluid have been derived using the Poisson bracket (PB) method 
by Son in Refs. 0, [ll- this approach one considers that the superfluid properties of a system arise from the 
spontaneous breaking of a continuous U{1) symmetry, with the appearance of a Goldstone mode. Since hydrodynamics 
is an effective field theory valid at long time and long length scales, the standard fluid variables should couple to the 
Goldstone mode. 

One defines the Hamiltonian of the system as a functional of the hydrodynamical variables and of the Goldstone 
field, and then postulates a set of PB among them to derive the hydrodynamical equations. The postulated brackets 
might be deduced from the canonical quantum commutators in the corresponding microscopic quantum theory. Some 
of them also follow from simple physical considerations. 

The equation of state of the system can be deduced from the microscopic theory and one finds that the differential 
equation for the total pressure of the system is given by 

dP = SodTo + nodfio + ^V^-d{d^ipf , (Al) 

where V is a parameter depending on the microscopic variables. As an example in complex scalar theories V is 
proportional to the expectation value of the scalar field 0, [1] ■ 

After a Legendre transformation of the pressure, one has that the energy density is given by p = sqTq + uq^o — P. 
Imposing Poisson brackets conditions on entropy density, particle density and momentum density one arrives at the 
hydrodynamical equations 0, [l] 

dpTP" = , TP" ^(p + P) uPu" - Pr/P" + V^dPipd^ip (A2) 

dp {nou" - V^dPf) = (A3) 

dpisoun = (A4) 

u^'d^.ip + /io = . (A5) 

In this formulation of the superfluid hydrodynamics there is a clearer interpretation of both the stress-energy tensor and 
current, as being due to the sum of the normal fluid part and the coherent motion of the condensate (the superfluid). 
Indeed both these quantities have contributions proportional to the hydrodynamical velocity u'*, associated to the 
normal fluid, and to the parameter V. Notice that the entropy has only contributions due to the normal component, 
as it should be. 

While in principle this formulation of the hydrodynamics looks very different from the one presented in Section |TT] 
they are equivalent. Indeed it is possible to relate the various variables of the two approaches by means of the following 
equations 



quiva 



A = p + V^{d^f , p^, = -dpip, (A6) 
= Sou" , nP = nouP - V^d^f , (A7) 

8^ = — [{TqSq + pono) Up + nodpip] . (A8) 
So 

Even if the two approaches are very similar they differ in one subtle point. While in the PB approach the value 
of uP-dp(p is constrained, in the convective approach this is considered as a free parameter. Instead one fixes the 
norm of the gradient of the field ip, requiring that at T = one has dpLpd'^Lp = p^. Imposing the two normalizations 
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simultaneously over-constraints the system, reducing the degrees of freedom that the problem should have. Let us 
also comment that the quantity n ^ /Iq. Indeed should be considered as the chemical potential associated to the 
normal component of the superfluid evaluated in the comoving frame. Instead /i is the norm of the vector fi^. 
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